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Abstract 

\^ ' We show that in certain hmits the (l+l)-diniensional massive Thirring model 

cn , at finite temperature T is equivalent to a one-dimensional Coulomb gas of charged 

particles at the same T. This equivalence is then used to explore the phase structure 
^^ ■ of the massive Thirring model. For strong coupling and T ^ m (the fermion mass), 

On , the system is shown to behave as a free gas of "molecules" (charge pairs in the 

Coulomb gas terminology) made of pairs of chiral condensates. This binding of 

chiral condensates is responsible for the restoration of chiral symmetry as T — > oo. 

In addition, when a fermion chemical potential ^ 7^ is included, the analogy with 
q"| a Coulomb gas still holds with fi playing the role of a purely imaginary external 

i-^ ' electric field. For small T and fx we find a typical massive Fermi gas behaviour 

for the fermion density, whereas for large fj, it shows chiral restoration by means of 
k> i a vanishing effective fermion mass. Some similarities with the chiral properties of 

5-H ' low-energy QCD at finite T and baryon chemical potential are discussed. 



1 Introduction 

The massive Thirring (MT) model in two space-time dimensions has been widely studied 
as a toy counterpart to low-energy QCD, since it does not include many of the complica- 
tions arising in 3-|-l dimensions. Amongst the features shared by the MT model and QCD 
is hadronisation (bosonisation). In this primitive version, the MT model is equivalent to 
the sine-Gordon (SG) model, both at zero and non-zero temperatures [||, ^, |^. Viewed 
as a non-linear signia model (NLSM) in 1-1-1 dimension for a single Goldstone boson with 
an explicit symmetry breaking term, the SG model mimics the chiral Lagrangian for 
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low-energy, strongly coupled QCD (whose role is played here by the MT model) where 
the lowest energy excitations of the vacuum are the Goldstone bosons. In addition, the 
solitonic excitations of the SG model (kinks, for brevity) correspond to Thirring funda- 
mental fermions and hence they are the analogue of QCD baryons or, alternatively, the 
NLSM skyrmions H- 

This paper is the second in a series concerning the statistical mechanical features 
of those models. In [^], the MT/SG system and their equivalences (bosonisation) were 
analysed at finite temperature T > and non-zero fermion number chemical potential 
/i 7^ 0. It was shown that when there is a non-zero net fermion number in the MT 
model, a topological term arises in the dual SG model which counts the number of 
kinks minus antikinks 0]. Physically this term represents the thermal bosonisation of 
fermions into kinks and it reflects the existence of pure fermion excitations in the thermal 
bath. The same term appears in the exactly solvable massless Thirring model ||^, and 
a similar contribution for a nonzero baryon chemical potential in the low-energy QCD 
chiral Lagrangian was obtained in p. However, the massless Thirring model, dual to the 
Schwinger model, has a much less rich and relevant structure if we have QCD in mind, 
and we shall not consider it here, except as a limiting case. 

The purpose of the present paper is to analyse the thermodynamics of the MT model, 
calculating physical quantities such as the pressure, the fermion condensate (the order 
parameter of the chiral symmetry) and the fermion number density. We believe this to 
be important for obtaining a better understanding of crucial physical phenomena such as 
the QCD chiral phase transition, not only at finite temperature but also at finite baryon 
density 0. 

We now introduce the model. In Euclidean space-time with metric (-I-, +) the La- 
grangian density of the MT model is 

Here ijj is a two-component fermionic field, the subscript denotes bare quantities and 
jfj,{x) = '4!{x)'y^'ip{x) where the Euclidean 7 matrices may be found in 0. For g"^ > 
forces are attractive and the theory is super-renormalisable. For —1/2 < g'^/n < 0, 
further renormalisations need to be carried out as discussed in reference P for example, 
whereas for g'^/ir < —1/2 the theory is no longer renormalisable [Q. Throughout this 
paper we take g^ > unless otherwise stated. With QCD in mind we are interested in 
strong coupling. In this we are aided by the general result that the natural parameter in 
which to express results is (1 + g'^/n)^^. This is understood by the explicit nature of the 
duality with the SG model, with Lagrangian density 



CsG[4>] = 7^d,(Pd'^(^-^cosX(P 



-d.(l)d^(t) - ^° 
2 ^ ^2 

where (/) is a real scalar field. Recall that these models are equivalent only in the weak 
sense, i.e. for vacuum expectation values and thermal averages, provided that their renor- 
malised constants are related through |I| 

(1.1) 
(1.2) 

where p is the renormalisation scale and m the renormalised mass at that scale. Thus 
large positive g'^ corresponds to small A^, for which approximations can be well controlled. 
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With QCD in mind, we will be interested in the chiral properties of the MT model as 
well as of physical quantities such as the pressure. The chiral U{1) transformations are 
ip —>■ exp[ia'y^)ijj for the fermions and — > + a/A for the bosons with a real arbitrary. 
The massless Thirring model is chiral invariant but the fermion mass term breaks that 
symmetry explicitly. Similarly the free boson theory is chiral invariant, whereas the 
cosA0 term in the SG model breaks it. Both models have still a residual Z invariance, 
corresponding to the choice a = 27m with n any integer. It is important to remark that 
in two space-time dimensions, the breaking of a continuous symmetry (f/(l) — > Z in 
this case) cannot be spontaneous [|T^]. The U{1) and Z symmetries are, respectively, the 
counterparts of the chiral and isospin symmetries for QCD, a^ and A playing the role 
of the pion mass squared and the inverse of the pion decay constant respectively in the 
effective chiral Lagrangian to lowest order, i.e. the NLSM |]11|. These generalities apart, 
the identification ( [L.l|) which, for strong coupling becomes Xg ~ 27r, is the only attribute 
of the SG theory that we shall use. Depending on which is most convenient, we shall 
switch between the use of g and A. This is not to say that the SG model has no interest 
in its own right. It fact, it is an ideal testbed for summation schemes for the pressure 



and, as such, will be considered elsewhere |jT2| . It is also of direct relevance to the study 
of kinks in Josephson junctions ||13 . 

In order to make progress we comment on a further property of the 2D MT model 
that also has a counterpart in QCD: at T = and /i = it is also equivalent to a 
2D classical statistical-mechanical system which consists of non-relativistic particles of 



charge ±g (a Coulomb gas) at a temperature Tcg H? |14|. Before clarifying the form of 
this analogy, we comment that one of the aims of this paper is to understand whether 
a similar equivalence holds when the MT model is heated to a non-zero temperature T 
and for /i 7^ 0. If it does, we will then be able to use the simpler Coulomb gas system to 
explore its phase space. It has not yet proved possible to perform similar calculations for 



QCD, as represented in |T3, in which the Coulomb gas is composed of monopoles. 



Indeed, one of our conclusions is that at high temperature T, the MT model can be 
equivalent to a one dimensional Coulomb gas of particles at the same temperature T. 
This is not the case for the T = MT model, for which 

J-CG 



A2 

leading to a Kosterlitz-Thouless (metal-insulator) transition at A^ = Svr [|14|. The ID 



Coulomb gas model has been solved exactly fT^, [T^ and we use those results to extract 



information on the behaviour of observables such as the pressure and chiral condensate 
of the MT model. Furthermore, when /i 7^ in the MT model, we show that the analogy 
with a ID Coulomb gas still holds with /i playing the role of a purely imaginary external 
electric field. 

This paper is organised as follows. In section ^ we analyse the conditions under which 
the equivalence which a ID Coulomb gas holds. Having specified this "dimensionally 
reduced" regime, the exact link between the parameters of the different models is then 
specified in section | where the pressure of the MT model is also calculated. Section |^ 
is concerned with the fermion condensate. We see that for large g and T ^ m, the MT 
model behaves as a free gas of "molecules" made of pairs of chiral condensates. This 
binding of chiral condensates is responsible for the restoration of chiral symmetry as 
T — » 00, and no phase transition takes place. The effects of a fermion chemical potential 
/i 7^ on the pressure and fermion density are discussed in section ^. 



2 Dimensional reduction of the T > MT model. 

At nonzero temperature T and zero chemical potential fi = the MT model partition 
function is given by [^ Q 
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(2.i: 



Here Zq{T) is the partition function for a two-dimensional free massless Fermi gas, which, 
ignoring irrelevant vacuum terms, is given by 



Zo{T) = exp 



ttLT 



and the function F2„(A, T,L) is 



F2n{X,T,L) = il / dTj I dq, n [Q\x, - Xk) 

j^l^O Jo i^^j 



ejekX'^/i-K 



(2.2) 



Herexj = {Tj,qj), (3 = T ^, 
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n + 1. 
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and L is the length of the system. We will eventually take the thermodynamic or L -^ oo 
limit in all our results. It must be clarified though that the L ^ oo limit will be taken 
always keeping (3 finite. In fact, we will see explicit examples below in which the L ^>- oo 
and T — > O"*" limit do not commute. Finally, the Q variable in ( |2.2|) also appears in the 
finite temperature free massless boson and fermion propagators P| and is given by 



Q^q, t) = smh(^^^— -) smh(- ^ ^ 



P 



P 



(2.3) 



betraying its conformal origins. As discussed in 0, the integrals in ( p. ID are convergent 
for g"^ > 0. 

Equation ( |2.1| ) will be our starting point here. First, notice that 



Q'iq,r) 
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(2.4) 



The key observation is therefore that if we were allowed to replace the Q^ functions in 
( p.l|) by their asymptotic values ( |2.4|) then, with appropriate definitions, (|2.1|) would 
resemble the grand canonical partition function of a one-dimensional classical gas of 
charged particles with positions labelled by g^. Remember that the two-particle Coulomb 
potential for charges ±a at points gi and g2 on the line is V^ oc ±cr^|gi — q2\- This ID 



Coulomb gas system was studied a long time ago [jT6|, |T^ and it can be solved exactly. 
We shall explore this analogy later on in section ^ and use the exact results to calculate 
the pressure and chiral condensate of the MT model. 

Before doing so, however, we need to establish the conditions under which the replace- 
ment given in ( p.4|) — which we will denote dimensional reduction (DR) — can be safely 
made. This DR regime is the 2D analogue of more complicated situations analysed in 



the thermal field theory literature, which roughly works for high temperatures and large 
distances [[l|, |l|, |20| . 

Clearly the regions in the integrand ( |2.1| ) where the replacement ( p.4| ) is not allowed 
are those where \qj — qk\ < P/t^- In those regions the dominant contributions to the 
integral come from Xj ^ x^ with e^efc = — 1 since 

Q\q, r) -^ {nT)\q' + r') for (g, r) ^ (0+, 0+). 



Thus, on the one hand, we expect that for high enough T the contribution (arising in 
the denominator) can be neglected. On the other hand, the above contributions become 
more important as g^ decreases, so that one may also think that for large enough g"^ 
the approach could be equally justified. Bearing these considerations in mind, we shall 
analyse the limits of high T and large g^ separately in sections 2?T and |2.2| . 

Before proceeding, two comments are in order regarding the expression (|2.1| ). First, 
every term in the n-sum picks up an overall j3L factor. In other words, the number of 
two-dimensional independent variables in the integral is actually 2n — 1. This is most 
easily seen by changing variables to 



Zi =Xi -X2,Z2 = X2 -X3, 



Z2n-1 — X2n-1 — X2n, Z2n — X2n 



(2.5) 



so that in (O 



fe-i 

Xj Xjs / ^ Zi 

i=j 
and therefore the integrand is independent of Z2n, which yields the jSL factor. Notice 
that this is typical of closed loops in perturbation theory [^. Order by order one has 
to consider all possible connected closed diagrams for the partition function and the PL 
factor is just the consequence of total energy-momentum conservation. One should bear 
in mind that the physically relevant object is not the partition function but the pressure, 
defined in the thermodynamic limit as 

P= \im —log ZMTiL,T) (2.6) 

which behaves as an intensive quantity. 

The second comment concerns the scale dependence. The partition function ( p.l|) is 
scale independent (and so is the pressure) since the explicit dependence on the renormal- 
isation scale p is exactly cancelled by the implicit dependence of the mass m(p) (see 
for details). Thus, unless otherwise stated and whenever dealing with scale- independent 
objects in the following, we will choose for convenience p = m, the renormalised mass of 
the Thirring fermion. 



2.1 High T limit 

Let us rescale qj -^ qjP and r,- -^ TjP Vj in (|2^ )Pi Then all the relevant T dependence 
is now outside the integrals in ( |2.2| ) since in the thermodynamic limit rj = L/jS —* oo the 
integrals are finite as long as A^ < 4tt — excluding of course the rj overall factor mentioned 
above. Then, choosing p = m, the series in ( ^.1| ) yields effectively a perturbative series 



in {m /T 
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^For much of this section it is more convenient to work with A rather than g in the first instance. 
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Figure 1: The function A/2 (A) 
with /2n(A) independent of 77 in the thermodynamic hmit and given by 
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/2n(A) = n / rf^? / dz] n 



i=k 



where Zj = (z°, z^) is defined in (|2.5|). It follows from ( p.3| ) that Q {Pq, (3t) is independent 



of /3. Therefore, to leading order in m/T, we can write the MT pressure as 



PMT(A,T) = P^;r(A,T) 



6 /m\^~^^/'^'^ 
l.U^) A/.(A) 






(2.7) 



where the asym superscript denotes the value obtained by replacing the Q's by their 
asymptotic values in (|2.4|), and A/2 (A) = /2(A) — /2*^'"(A) with 



/2(A) 



dq / rfr Q\pq,pT) 
JO L 



(2.8) 



and 



/2^'"^(A) 



2A2/27r+l 

A2 



We have evaluated numerically A/2 (A) and the result is plotted in Figure Q. We see 
that it remains of C(l) until very close to the limiting case A^/47r<l, or g^>0, where it 
diverges. Thus, the relative error for T ^ ?ti is at least of order (m/T)'^ and we can 
therefore conclude that for T ^ m, DR is valid for g'^ > 0. 



2.2 Strong coupling (large g'^) limit. 

We would now like to study the pressure for g^/n ^ 1 (A^/47r ^ 1). The first observation 
is that the integrals in ( |2.2| ) diverge if we set A = and then take the limit L ^ 00. 



^Notice that this is a scale dependent condition 



Therefore we shall keep L finite and then take it to infinity (keeping (3 finite) only when 
the pressure (|2l^ ) is calculated. Then, expanding (|2.2|) in \^ /An one obtains 



nl3X^ p J /■! 
7tL Jo Jo 



F2n{\,T,L) = {(3Lf^ I 1 - '-^ I'dq Idrlog Q\[3q,[3T) 
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F2.(A, T, L) = F.T'^iX, T, L) - n{PLf--I{r^) + O 



An 



An 
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dq dr log Q^{pq,pT) +2r]\og2-nr] 






(2.9) 



X' 



2n 
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where 



Zt^y-{T,X,L) - ^nv)Z^iT)-I'oiz) + 0^^j 






and /o(-2) is the modified Bessel function of zeroth order. On the other hand 

oo . 1 

Zmt(T,A,L) = Zo^(T)5: - 



n=0 



l\Vz^2n 



n! 



27 +^(£)-^°^^^)^°(^) + ^(£^ 



(2.10) 



(2.11) 



and Zl^j^{T, A, L) has the same leading order in X^/An as above. Therefore, we can write 
for the pressure 



PMT(T,A) = P^;r(r,A)-lim — 



Now, for large z |22 



e^ /I 

y'2nz \z 



8w ™ /o(^) 



for z^ 1 



O 



An 



so that taking L — >• oo we finally get 



Pmt{T, A) = Prr(T, A) - ^/(ry - oo) + O [^) . 



(2.12) 



Numerical analysis of the function I{ri) shows that it clearly vanishes as 77 — i> oo. There- 
fore we have APmt = O(A^). Although it might seem that this result is valid for any 
T, one has to be extremely careful when approaching the limit T -^ 0~^. In fact, we 
see that expanding ( p.lO|) in A yields logarithmic factors logT/m hidden in the (9(A'^) 
in ( p.l2|) , but which will show up at next-to- leading order. Therefore the above results 



■^Note that, in turn, from ( ^.11 ) we find the asymptotic Umit Pmt — fT^ + rn^ for small A, which we 
will recover in the Coulomb gas approximation in section ^ 



should not be trusted for temperatures T <^ m. This is just a consequence of the non- 
uniformity of the T —>■ 0'^ hmit or, in other words, that the L ^ oo and T ^ 0+ do not 
commute in general Hence we will consider that our approach is justified for large g^/n 
and T>m,. Recall that the presence of the log T/m factors in the A expansion may be 
also troublesome for very large T, which gives a hint that the high T limit may not be 
compatible with perturbation theory, as has been noticed in more complicated situations 
In our case though, we have a well defined high T expansion, namely our expansion 



in (m'^/T'^Y^^ ^^^ discussed in section p.l| , which we will indeed be able to resum in the 
DR regime. 

To summarise, the results of this section are that the range of validity of the DR 
approach which we will be considering for the pressure is 

T:^m , / > 0, 

T>m , ^>1. (2.13) 

TT 

We now turn to the ID Coulomb gas system and specify exactly the link between it and 
the MT model. 



3 The MT partition function as a Coulomb gas 

Consider the partition function for a one-dimensional neutral classical non-relativistic gas 
of charged particles. In particular, take A^ positively charged particles and A^ negative 
ones, where the magnitude of the charge will be denoted by a. In one dimension, these 
charges interact via a Coulomb potential proportional to ±o"^ and to the distance between 
them on the line. This system can be also interpreted as uniformly charged plane sheets 
moving along the direction normal to their planes |T6|, IT^. Then, the partition function 




n{z, e, a,L) = J2 j^ I n L "^li 1 ^^p 



i<i<«<2iv 



in the grand canonical ensemble at fixed temperature 6, fugacity z and length L is HTq , [T7|| 

(3.1) 

where e, = 1 for i < N and e^ = — 1 for i > N. The fugacity is related to the chemical 
potential fi hj z = V^ttMO exp 9~^n, where M is the mass of the particles[|. However, 
we will keep z instead of /i here, so as not to confuse this chemical potential with the 
fermion chemical potential we will introduce in section ^. In the thermodynamic limit 
the pressure and the mean particle density are given by 

Piz,9,a) = lim -\ogn{z,9,a, L) 

n{z,e,a) = = 6 z—P{z,e,a) 

where <C ■ ^ denotes the statistical average in the above ensemble. We therefore realise 
that by making the replacement ( |2.4| ) in ( p.l| ) and setting p = m, one can writeQ 

Zmt{T,L) = Z^{T,L)n {z = — i^—j ,^ = T,a = -=,L . (3.2) 



"'This definition of fugacity differs from tfie usual one by tfie factor of \/M. We have absorbed this 
factor into the definition here as M plays no role in future considerations. 
'"^ Again we work first with A and translate to g later. 
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This is the central equivalence of this paper. Its utility lies in the fact that the 
above classical problem was solved analytically in references 



Mm 



There it was shown 

that for small 6 the ID Coulomb gas behaves as a gas of free "molecules" , made up of 
H — charges pairs bound together (since the mean kinetic energy is much smaller than 
the mean potential energy). The main features of this phase are: i) The pressure is 
small compared to 27ra'^, which represents the pressure between H — charges, and ii) The 
probability density of finding a H — pair within a given distance r is much bigger than 
that of finding a ++ pair for distances r less than a typical "molecule" size, which on the 
other hand is much smaller than the typical inter-particle distance. On the contrary, for 
large 6 (when the mean kinetic energy is much larger than the mean potential energy) the 
charges are completely deconfined, forming an electrically neutral "plasma" of 2N free 
particles, where the pressure is higher than in the "molecule" phase. The crucial point in 
our case is that, from (|3.2| ) we see that the unit charge grows with temperature (a oc T) 
and therefore the above phases are going to be reversed (high T "molecule" phase and 
low T "plasma" phase), as we shall now see in detail. 
Let us first quote the result for the pressure in p!7| p| 



P{z,9,a) = 2na'^-fo 



zO 



27ra^ 



where 7o(-2) is the highest eigenvalue of Mathieu's differential equation 



+ 2z cos ( 



y{.<P) = iy{.<P) 



(3.3) 



withj/(0+27r) = y(0). Notice that (I 



differs slightly from the notation customarily used 
for Mathieu's equation (see for instance p2|), which is recovered by setting 7 = —a/4, 
z = — g/4 and (p = 2v. Therefore, from (|3.2|) we obtain for the full pressure of the MT 
model in the DR regime: 

._, (3,4) 



Pmt{T) = t:T- + Pc(T, A) 





with the "Coulomb gas" pressure being given by 



Pc{T, A) 



-70 



m 



A2T2 \mj 



A2/47r' 



Or, in terms of g"^, 



Pc{T,g) 



27rT2 



l + ^Vvr 



7o 



_rri^( g^ 



2T\ (i+sV'r)-'' 
m J 



(3.5) 



The eigenvalues and eigenfunctions of the characteristic problem ( |3.3| ) are well known 



and tabulated. We refer to [^] for their general properties and quote here only those 
relevant for our purposes. For instance, the asymptotic limits of 7o(-2) for very large and 
very small z are 



7o(l) ~ 2z^ + 0{z^) for z < 1 
7o(z) ~ 2i-A/i+0(l) for z > 1. 



(3.6) 
(3.7) 



®Note that in ||l^, the pressure and other quantities are calculated in units of 9/211 = 1 with a integer. 
Changing variables in (3J) as qi — > qiO/iira'^, we have Vt{z,0,a,L) = 51(z0/27r(T^, 27r, 1, 27r(T^6'~^L) so 



that Edwards and Lenard's results can easily be rescaled for our purposes. 
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The relevant parameter setting the quahtative behaviour of the system is the argument 
of 7o in (U), 

1 + ^ — ■■ (3-8) 



AttT^ \ TT J \ m J 

for small z, the system is in the "molecule" phase, whereas for large z it is in the "plasma" 
phase. Let us then consider three separate hmiting cases: 

i) 5'^ ^ ^ ^"^^ T ^ mg/2'K\ This is the very high T regime and we have i ^ 1, 
so that equations (^]6|) and (|3.5| ) give Pc ~ w? {w? g^ / Ati'^T'^) , i.e. the Coulomb 



pressure vanishes at high T for large g . This is the "molecule" phase described 
in 



16|| in which the Coulomb charges tend to pair, thus lowering the pressure. As 
noted above, this "unexpected" behaviour of high temperature binding is just a 
consequence of the fact that our coupling "charge" a in (|3.2|) increases with T. 
Besides, when adding the free term in ( p.4|) , the total pressure is an increasing 
function of T, behaving in the high T (and large g"^) regime as a free gas of massless 
fermions. This is indeed the first indication that the "molecule" phase corresponds 
in fact to a phase of asymptotic chiral symmetry restoration. The reason is that 
we expect that if the chiral symmetry is restored, the system behaves roughly as 
the massless case which is chiral invariant. In other words, we would expect the 
effective fermion mass to vanish in the T —>■ oo limit when the system tends to 
restore the chiral symmetry. 



g"^ ^ Ti but T <C mg/2T:: Now equation ( |3.7| ) gives Pq ^ m? — nmT/g. Since 
i 3> 1, now we are in the "plasma" phase. Notice that this same result, to leading 



order, could have been obtained through our analysis for large g in section 2.2 



c.f. equation (|2.11|) . Again we note that to (9(A^) in the expansion of Pc there 



is a contribution proportional to X'^m'^ log{T/m) which prevents us from going to 
T <^ m. This is another reason why in the small A (large g) regime we can only 
trust DR for T>m, which is the region covered in this second case. 

iii) < g^ '^ 1 and 27rT 3> m: we have now z ^ m/2TiT . Again, z <C 1 and we are 
in the "molecular" phase. Remember that in this region we can only trust DR for 
T ^ m. Thus Pc ~ m^/n so that the pressure tends to a constant value for high 
T. Actually, if we were able to push our results formally to g^ < (A^ > in), 
the pressure would start growing for large T. Clearly, our previous calculations 
in the MT model are not justified in this case, which at least would require extra 
renormalisations. However, from the viewpoint of the Coulomb gas only, there is 
no particular problem with taking g"^ negative. In other words, the dimensionally 
reduced theory is UV finite. Notice also that, as we make g"^ more negative, the 
Coulomb correction to the free boson gas term becomes more and more important. 
In particular, if we take 5'^/vr< — 1/2 f\, Pc would grow quadratically with T, thus 
being of the same order as the free contribution. 

In Figure ^ we have plotted Pc{T) in units of m for different values of g^. Notice that, 
as we have said before, we cannot claim that curve iv) corresponds to the MT model in 
the DR regime. We are simply extrapolating the Coulomb gas results to that point. As 
for curve iii), we should intend it as the limit of the dimensionally reduced MT model for 
very small but positive g^, where we can still trust our DR approach, as long as we look 

^Recall that this corresponds to A^>87r and thus to the Kosterhtz-Thouless transition in the zero- 
temperature theory. 
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Figure 2: The Coulomb pressure Pc(T,g). 

only to the T/m ^ 1 tail of that curve. Note also that Pc{T) is always a continuous 
function of the temperature and thus we do not see a phase transition in T, consistent 
with being in two dimensions. In Figure ^ we plot the Coulomb pressure and the total 
pressure for g'^/ir = 1 in order to estimate the size of the "Coulomb" gas correction to 
the free gas. 



4 The chiral charges and "molecules" 

In the previous section we have exploited the analogy of our dimensionally reduced MT 
partition function with that of a Coulomb gas on the line in order to calculate the pressure. 
In this section we will show that this analogy can be extended further, i.e. that it also 
works in the same way for other observables. In particular, we will concentrate on thermal 
averages of products of the operators 

since these account for the chiral properties of the system. Under the chiral transfor- 
mations discussed in the introduction, a±{x) -^ exp(±2m)cr±(x); in other words, the a± 
operators have well defined ± chiral charge. Below we will show that in our Coulomb gas 
analogy they play the role of the charges, and the chiral invariant combination a^a^ will 
then represent a "molecule". Thus by forming "molecules", the system tends to restore 
the chiral symmetry, i.e. it behaves as the massless theory where only o"_|_o"_ combinations 
are allowed. 

We will show below that the a± correlators can be related with the ID Coulomb gas re- 
duced density functions discussed in |]1^. These functions are defined as follows: f±{x)dx 
denotes the probability of finding a ± charge in the element dx; /+_(xi, X2)dxidx2 is the 
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Figure 3: The total pressure and the Coulomb pressure for g'^/ir 



joint probability of finding a + charge in dxi and a — one in dx2, and so on. It turns out 
that such density functions can also be exactly calculated Ijl^. We will only consider here 
the one-point and two-point functions, which will provide information on the relationship 
between Coulomb "molecule" pairing and the chiral symmetry, as explained above. 



4.1 The fermion condensate and chiral symmetry vs "molecule" 
pairing 

The first correlator we will analyse is just the one-point function, i.e. the condensate of 
a±. Observe that ^ o"+ ^ = ^ o"_ ^ since the MT model is invariant under ■?/'+ <-^ ip^ 
where 'iIj± = (1/2) (1 ± 'j^)ip are the right and left-handed projections of the spinor field 
{(j) — >■ —(J) in the SG model). This is a parity transformation and has nothing to do with 
the chiral transformations we have been discussing here. We shall implicitly make use of 
this -|- ^^ — symmetry in the correlators throughout this section. 

The scale invariant fermion condensate m -C ipip ^ is the order parameter of the 
chiral symmetry, exactly as the quark condensate in QCD. An immediate consequence 
of the absence of a phase transition in 2D is that the fermion condensate cannot vanish 
strictly at any temperature. On the other hand, for high temperatures T ^ m we would 
expect that the mass scale m becomes irrelevant and thus the system would restore 
the chiral symmetry, as indeed happens in the QCD chiral phase transition. This was 
already suggested by our previous analysis of the pressure. We may therefore expect the 
condensate to become smaller at large T but never to reach zero. However, we do not 
expect chiral restoration for all positive values oi g^. An indication that this may actually 
be the case is the following. In the g^ — > 0+ limit, the system should behave as a free 
massive fermion theory. However, it is not difficult to see that for g = the condensate 
behaves for T ^ m, as ^ ipip ^^ (m/vr) log(T^/m^), so that there is clearly no chiral 
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restoration in that case. Therefore we expect that chiral symmetry restoration is lost as 
the value oi g^ is decreased. Indeed, we have seen that behaviour already for the pressure, 
where the Coulomb correction to the massless term did not vanish for large T and for 
very small g"^. 

We will now show that in the DR regime the fermion condensate can be calculated 
exactly, thanks again to the analogy with the Coulomb gas. In fact we do not need to 
appeal to the density functions in this case since in the thermodynamic limit L ^ oo the 
system is translationally invariant so that 



- d 

m <C ipip > = m- — Pmt{T). 
am 



(4.1) 



Notice in particular that this implies that the regime in which we are allowed to replace 
the Q's by their DR values for the condensate will be the same as for the pressure, i.e. 

Thus we only need to differentiate once in the expression obtained for the MT parti- 
tion function in the DR regime to get the condensate. However, in order to clarify the 
procedure we will follow for the two point correlator and to understand better the anal- 
ogy between Coulomb charges and chiral operators, let us relate the condensate with the 
one-charge density functions f± analysed in [l7|. Those functions are given by /+ = /_ 
with 



U{z,9,a,L;X) 



]_ oo ^2N /27V 



E 



^ ^^l {N\)^ 



N= 



X exp 



2N L \ N 

n/ dqAY.Siq,~X) 



l<i<«<2Af 

^2N 



Qj\ 



n 

X exp 



E 



nJiNKN- 



^(nf^..)%.-^) 



271(7^9 ^ Yl (^i(^j\<li -?jl 

l<i<«<2Af 



(4.2) 

where Q is the partition function of ( |3.1| ). We see that in the L — ;> oo limit, a simple shift 
Qi ^ Qi + X ^i = ^, ■ ■ ■ , ^N ensures that /+ is independent of X. In turn, notice that 

U{z, 9, (r,L) = j £ dXUiz, 9, a, L; X) = ^z^ logQ^z, 9, a, L). (4.3) 

Now compare ( [4.2| ) with the expression obtained in for the condensate: 



m -C ipip S> 



Zmt(T) 

n+l ^ 



n=0 
■L 



nUn + 1)! 



pm (T 



n / d^^ i d<ij n [Q"i^^ 
j=i ■^'j -^^ k<j 



•^k/ 



ejefeA2/47r 



Thus, shifting n ^ n — 1 in the above equation, replacing the Q's by their asymptotic 
values (|2.4| ), fixing p = m and comparing with /+ we obtain in the DR regime 



m < V-V^ > = 2T/+ z 



m 



2 ,!yTi.\^lA-K 



2T \mj 
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9 = T,a 



XT 



(4.4) 



Relation ( [4 .41 ) is very interesting indeed since it supports the idea that the role of the 
Coulomb charges is played here by the a± chiral correlators (in the DR regime). We will 
elaborate further on this issue below. 

As it was said above, the fermion condensate can be obtained either using the analogy 
with the f± functions or by differentiating the pressure. We should then be able to rewrite 
( [4.4| ) as ( [4.1| ), which would be a consistency check of our results. For that purpose, all 
we need to use is ( [4.3|) and the equivalence between Vt and Zmt found in the previous 
section. In doing so it is important to leave the scale p unfixed, since differentiating with 
respect to m is a scale-dependent operation. After multiplying by m the result becomes 
scale independent, and we can safely fix p = tti. Once that consistency check has been 
performed, let us write the final expression for the fermion condensate as 



m -C ipip 3>= 



m2 /2T\(i+^'/-)"' 
2 \m) 



7o 



m 



47rT2 



'-)(?L\ 



n 



m J 



(4.5) 



where 7o(-2) = d'~fQ{z)/dz and we have converted back to g^. Remember that unlike the 
partition function, there is no contribution to the condensate coming from the free Bose 
part since the free theory is chirally symmetric. Therefore by taking into account the 
asymptotic behaviour ( |3.6| )-( |3r^ ), we now analyse the same limiting cases considered for 
the pressure in the previous section: 

i) q'^/ti ^ 1 and T ^ mg /27c: we have m -C ipip ^ ~ {m? /2Ti'^){mg/TY. The 
condensate vanishes asymptotically and hence the chiral symmetry is restored as 

T — > oo. 

ii) g'^/n > 1 and T < mg /2tx: Now m < ^V^ > ~ m2(2T/m)(^+^'/'')"\ This is the 
behaviour for temperatures T>m. Notice again that if we insisted on extrapolating 
this behaviour to T = we would find <^ ^^ S> (T = 0) = 0, which is clearly 
incorrect since we know that at T = the chiral symmetry is broken. Once more, 
we remark that we should trust our result only for T>m due to the presence of 
logarithms in the A expansion. 

iii) < g^ <^ 1 and m <^ 2ttT: m -C ipip ^ ~ 2m^ j-n^ so that the condensate tends 
to a constant value at high temperatures (just as the pressure did), and there is 
no chiral symmetry restoration. As we anticipated before, we do not expect chiral 
restoration to take place for small g^ . In turn, notice that the exact g^ = ^ limit 
predicts that the condensate grows logarithmically for large T, while we get here 
that it remains constant for g^ small. This is another sign of the special character 
of the (7^ = case, where, as commented several times before, we cannot apply our 
DR arguments because of the extra UV infinities. 

In Figure |^, we have plotted the condensate as a function of the temperature for different 
values of g^ and we can see the different asymptotic limits discussed above. It must be 
stressed again that the g^ = case plotted in that figure should be taken as g^ very small 
but positive. 

From the above discussion the following picture emerges: for large values of g^, the 
system at high temperatures behaves like a gas of neutral "molecules" made of pairs of 
chiral charges (chiral neutrality) and the chiral symmetry tends to be restored continu- 
ously (i.e, no phase transition) and asymptotically as T increases. As we decrease g^, the 
"molecule" phase tends to disappear in favour of the "plasma" phase — as was already 
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Figure 4: The ferinion condensate as a function of temperature in units of m, for different 
values oi g'^ 

seen for the pressure. Here the a^ and cr_ correlators can be different from zero in chiral 
non-invariant combinations and thus the chiral symmetry remains broken even for large 
temperature. 

Consequently if we look now at the two-point correlators, we should see a tendency 
of the system to increase the cr+cr_ correlator against the o"+cr+ one in the "molecule" 
phase, i.e, for high temperatures and large enough g'^. That will be the purpose of the 
next section. 



4.2 Two-charge correlators and the screening length 

Let us begin by recalling the definition of the two-charge density functions in the Coulomb 
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(4.6) 



X exp 
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! (n // dq,) 6{q^^, - X)6{qN - Y) 



X exp 






(4.7) 



As in the single charge case, the above functions depend only on X — F in the L ^ oo 

limit (translation invariance). Also, notice that /_+ = /+_ and / = /++ while for 

brevity we have omitted the dependences on z, 6, a, L. Our task now will be to compare 
( [4.6D -( ^!7|) with the MT correlators at finite temperature; <^ (T+cr_ ^ and <^ cr+cr+ ^ 
respectively. Let us start with the H — correlator. In P], the following result was obtained: 



<Tc(7+(x)(7_(r)> 



Z^{T) /TV /'" /p\2 - /l\2 
Zmt{T) [pJ [2) h\n\) 
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^Jl Jo 'Jo 
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X^/Att' 
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where e^ = + for j = 1, . . . , n + 1, e^ = — for j = n + 2, . . . , 2n + 2 and Tc means contour 
ordering along C = [0, —i(3\. For convenience we have retained A in the exponents rather 
than g. Again, the above correlator depends only on X — F in the thermodynamic 
limit. Notice that we are using the same notation X and Y both for one-dimensional and 
two-dimensional variables; the meaning should become clear from the context. 

As for the a+a^ correlator, it can be calculated through the same procedure followed 
in for the H — correlators using the generating functional technique. Notice that this 
correlator is not chiral invariant, unlike the a+a- one, which makes its analysis slightly 
more involved, as explained in P]. Despite that, we obtain after renormalisation 



< Tca+{X)a+{Y) > 
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(4.9) 



where e^ = + for j = 1, 
Note that the ratio 



k<j 
. ,n + 2, e^ = 

i?(A,T;X) 



- for j = n + 3, . . . ,2n + 4. 

<^Tccr+{X)a_{0) > 
<Tca+(X)a+(0)> 



is scale invariant — this is the observable in which we are interested. However, before 
proceeding further, we need to clarify how DR works for the above correlators. We begin 

16 



by noting that DR is not allowed for all values oi X — Y in ( ^■8| )-( ^l9| ). In fact, performing 
the 5^'^^ integrations, a factor [Q^(X — y)]^^ ^^^ comes outside the integrals and it is clear 
that the replacement (|2.4|) will be allowed for that factor only for q ^ (3/27T, where q is 
the spatial component oi X — Y. Besides, our analysis for the partition function does 
not imply that even inside the integrals in ( [4.8|) -( ^l9D such replacement can be done. The 
detailed analysis (which follows a similar line to that of sections p.l| and p.2| ) can be found 
in the appendix. Here we simply summarise the result: the structure of the correlators 
( [4.8| )-( ^^ makes it possible to reduce them dimensionally (inside the integrals) in the 
large g^ (small A) limit but not in the large T one. If, in addition one works at large 
distances q ^ P/^n, the the Q's outside the integrals are also dimensionally reduced. 
Bearing this in mind, we shall proceed to evaluate these correlators in the dimensionally 
reduced regime by exploiting again the analogy with the Coulomb gas. 

4.2.1 Two-point correlators and the Coulomb gas 

Given the conditions for DR to work for the correlators, let us now compare ( [4. 81 ) with 
( [4.6|) and ( [4. 91) with (|4.7| ). For that purpose we replace the Q's inside the integrals in ( [4.8|) 
by their asymptotic values ( ^.4| ). In doing so, we have to remember that such replacement 
is not allowed in general for the term j = n + 1, k = 2{n + 1) (see the above comment). 
Next, shift n — > n — 1 in the sum so that, comparing with (|4.(j|), we get 

«rca+(X)a_(0)» = - — 

\AJ rn^[Q^{q,T)] ' 



m- /2r\^/^" „ „ AT 



i2 



^ ^-l^=2flmj •'-^•"-JS-'" 



where X = (g, r), we have set F = and we have also made use of (|3.2|) . It is important 
to point out that the overall factor T"^ /m? depends on p implicitly through m but does 
not have explicit p dependence (we cannot fix p for that factor), whereas the rest — i.e. 
the function /+_ — is scale independent and we have fixed p = m only in that part. 
Following the same steps for the ++ correlator we now find 



<Tc(T+(X)a+(0)> = 4 — -^^^^ 



A2/47r 



^2 ^2T\^ /^" „ _ AT 
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and therefore the scale-independent ratio of the two yields 

AV47r 



fl(A,r;A') = 2-^/-— — -^ ) :^-^, ^-4. (4.10) 



ma,r)f''' uJ.= '4,{^f"\e = T,a 
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Once more, the main advantage of comparing with the Coulomb gas is that the 
functions /++ and /+_ can be calculated exactly. They can also be related with the 



Mathieu equation (p.3|) as []17 
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Figure 5: The two-point correlator ratio R{T, g; q) for different values of T and g^. The 
numerical parameters for the above curves are: i)z = 0.225, (27rT)~^/go = 0.044, ii) 



z = 0.082, (27rT)-Vgo 
(27rr)-Vgo = 5 X lO'^ 



1.3 X 10-^ iii) z = 0.01, {2nT)-^/qo = 2.3 x 10-^ iv) 



0.005, 
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the L ^ oo limit has been taken, z = 6z/2'n-o'^, q = 2'Ka'^q/9 and ym{4>), 7m are. 



respectively, the eigenfunctions and eigenvalues of ( |3.3|) , with 7m-i > 7i 



m 



1,2, 



First, note that since 7„ 
above sum in the limit q - 



< for m > 1 and 70 > 0, only the Bq term survives in the 



00. Then, since 2/0 (0) 



yo[ 



we have /+± -^ (Bq) 



m 



Since in that limit we can also dimensionally reduce the 



that limit. In fact Bq = Z'Jq^z) 
Q^ factors in (|4.10|) , we get 

R{X, T; g -^ 00, r) = 1. 

After conversion to g, our results for R{g,T]X) are plotted in Figure ^ for different 
values of T and g"^. In all those curves, we have considered values q ^ (27rT)~\ so that 
we can neglect the r dependence and thus ( [4.10D becomes just the ratio f+~{q)/ f++{q)- 
We also recall that the /++ function can take negative values for smaller q, but we have 
only displayed in Figure |^ the region where R{T, g; q) is positive. 

Let us define the parameter go to be the mean inter-particle distance so that go = 

1//+, i.e. the inverse of the density of charges, regardless of whether they are positive or 

negative. Thus, 

2T 
go(A,T) = 

m <^ ipip ^ 

with m -C ipip ^ ill (EiH)- 

As in our previous analysis of the pressure and the fermion condensate, the relevant 
parameter here is z in ( |3.8| ). For very small z the system is in the condensed regime. 
Since we need to work at large g'^/ir for DR to work, small z means large T. In that 
regime, the H — correlator is much larger than the ++ one at short distances (but yet 
large distances compared to (27rT)~^), as a result of the tendency of the system to form 
cr+(T_ pairs ("molecules"). For large distances, R tends to 1. In the "normal" phase, 
i.e, when no "molecules" are formed, we expect that R starts approaching its asymptotic 
value near g = go and it does not grow very much below go- Conversely in the "molecule" 
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phase, R grows for small q and it approaches 1 much faster, defining a certain qcond -C go, 
which we interpret as the screening length or "molecule" size. For the curves plotted in 
Fig.|, we get i) qcond ^ l-3go, h) q^ond ^ 0.3go, iii) qcond ^ O.Olgo, iv) qcond ^ 0.0025go- 
One can clearly observe that the screening length decreases with the temperature. 



5 Dimensional reduction of the Thirring model at 
nonzero fermion chemical potential 

In 0, the issue of thermal bosonisation in the MT/SG system was also studied for /i 7^ 0, 
where /x is the fermion chemical potential associated with the conservation of the fermion 
number density tp'^oip in the Thirring model (the conserved quantity is the number of 
fermions minus antifermions). It has been shown that the SG model Lagrangian acquires 
an extra topological term, interpreted as /x times the number of kinks minus antikinks. 
The purpose of this section is to show that in the DR regime, the partition function of 
the Thirring model at T > and /x 7^ can also be related with the ID Coulomb gas. 

We start by recalling the result (in perturbation theory in m) obtained in for the 
partition function of the MT model at T > and /i 7^ 0: 

Zmt{T, /i; L) = Zm=o{T, fi; L)Zc{T, fi; L) 

where Zm=o{T, fi', L) is the massless Thirring result derived in [^ and is given by 
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(5.1) 



(5.2) 



and, as in ( p.2| ), e^ = +1 for j < n, ej = —1 for n < j < 2n. Notice that the above 
partition function is real, as it should, because of the symmetry of the integrand under 
the simultaneous relabelling gi <-^ qn+i, ■ ■ ■, In ^^ 'l2n- 

Once more, we need to analyse the regime in which we can dimensionally reduce the 
Q variables inside the above integral through ( p.4| ). It turns out that the analysis in 
this case is much simpler than in the case of the two-point correlators and indeed the 
conclusion is that DR works in the same regime as for the ^ = pressure, that is, ( [^.13 ). 
We quickly sketch the argument here: first, for small A one readily notes that the leading 
order is /i-independent since J2j ^j = 0. Second, for large T, one arrives to an expression 
identical to (|2.7| ), with /2(A) in ( |2.8| ) replaced now by 



/2(A,/i) 
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whose numerical analysis shows a similar behaviour than /2(A), i.e, its correspond- 
ing A/2(A,/i) remains of C(l) for A^/47r<l for different values of /?//. That is eas- 
ily understandable by realising that /2(A,/i) < /2(A) and /2*^™(A,/i) > 0, so that 
A/2(A,/i) < /2(A), which near A^<47r (where the error is bigger) is /(A^<47r) = 0{1). 
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Given therefore that one can work in DR also for /i 7^ 0, the Q's in ( |5.2| ) may be 



replaced by ( p.4|) . We then reahse that Zc(T,fi) in ( |5.1|) is now the grand canonical 
classical partition function of the Coulomb gas, the extra /^-dependent term having exactly 
the form of a purely imaginary external electric field term added to the Hamiltonian of 
the system. In fact, note that if we wanted to describe (A^ positive and N negative) 
charges on the line not only interacting among themselves with a Coulomb force, but 
also with an external electric field S, we should add to the Coulomb Hamiltonian the 
term 

2N 

n^n-a£j2eiqi. (5.3) 



i=l 



However, this is precisely the form of the /i-dependent term in ( |5.2|) , so that once the Q's 
are dimensionally reduced one obtains 

Zc{T, ^;L) = n\z = —-[ — ^e = T,a = -y=, S = ^; L 

2T \m J V47r Vvr 



with Q{z,6,a,S,L) the grand canonical partition function of the classical Coulomb gas 
interacting with an external electric field S through ( ^.3|) . 

The exact calculation of the pressure P{z,9,a,S) = limL-^oo^ogQ/ /3L was also per- 
formed in ||l^.0 In this case though, the result cannot be expressed in a simple way in 
terms of Mathieu functions as for /i = 0. Now it is given by 

PciT,f^) ^ hm -^\ogZciT,fi;L) = ^^^^(z) (5.4) 

L^oo pL i + g /IT 



where z is given in (|3.8| ) and the function 7(z) is defined implicitly as follows. Let 
G{z, 7, 7]) be the following continued fraction: 



7+4J7+4- 
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where 



S . /i 
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and define 



47ra 27rT' 

Q{z, 7, ^) = - [G{z, 7, V) + G{z, 7, -7])] . 

7 

Now let zo{'y,ri) be the smallest positive solution to the equation 

0(^0(7, ^), 7, ^) = 1- (5.6) 

Then 7(5) is defined implicitly by identifying z"^ = zqI^jT]). 

First note that for /i = the above result reduces to the solution in terms of Mathieu 
functions we have found in section |^. In that case equation ( |5.6|) defines precisely the 
condition that the two parameters z and 7 of Mathieu's differential equation (|3.3|) should 



®Note that the ensemble used in llq] was the canonical one at fixed N and P (so that L is varying, for 
instance by means of a moving piston) rather than the grand canonical ensemble at fixed z and L being 
used here — the former are the conjugate variables of the latter. However, in the thermodynamic limit 
the results obtained in the different ensembles should be equivalent, after the standard identifications. 
Namely, for N — > 00, L ^ 00 with N/L fixed, we can identify N with N and L with L, where the bar 
denotes averages in the corresponding ensemble. 
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satisfy in order that the solution is 2tt periodic. Therefore for /i = 0, we have just 



7 = 7o. In [^ it was shown that equation (|5.6| ) has always a solution. It is worth 
mentioning that in our case, due to the purely imaginary character of t], the convergence 
of the different integrals analysed in |jl6| and hence of the continued fraction ( ^.5|) is 
automatically ensured. This should be contrasted with the real rj case where 7 and rj 
have to satisfy certain constraints to ensure convergence. In addition, notice that for 7 
real, Q in ( |5.5| ) is also real, which means that the pressure is real, as it should. 

Equation ( |5.6|) can be solved exactly in the asymptotic limit 5^1 and 7^1, where 
we know there is condensation into "molecules". For 2;,7 -C 1, Q{z,'y,ri) in ( [5. 51 ) just 
becomes 



7(1 - At]"^) 
so that 

- irj^ (5.7) 



z'^ ^ 



7UJ 



2 

is the smallest solution to Q = 1 and defines 7(i) in this limit. Therefore, using ( p.4| ), 
replacing the values for z and 77, and eliminating A^ for g"^, we find in the ;z ^ 1 limit 
(condensation regime) 

4n i^2/^\ .omx 2(l+gV^)-' 1 



P,(r,,).!!!%±^n 



47rT2 \mj 1 + /j'^/tt^T^ ' 

One can check that for fi = the above result reduces to the one we have found in 
section ^ in that regime. Notice that the pressure not only vanishes for large T {T ^ m) 
but also for large /x at large T {ft ^ T ^ m). That is, chiral restoration also takes place 
for large fi (in the large T limit here). This is confirmed by differentiating the above with 
respect to m to get the fermion condensate, which clearly vanishes as fi —>■ cxo. On the 
other hand, differentiating the pressure with respect to fi yields the net averaged fermion 
density, 

piT,fi) = hm J « V^7V » = ^PiT,fi). 

L^oo L Oil 

Thus, we get for the total fermion density in the condensation regime 

where the first contribution is the massless result obtained in |^. The above result for the 
fermion density deserves some comments. One notes that in ( p.8| ) the Coulomb correction 
to the massless case is always negative and very small in this regime, since the above 
result is valid only when T ^ m. Therefore as T ^ 00, only the massless contribution 
remains since the first term is T-independent. This is nothing but another signal of chiral 
restoration for large T: in that regime there remain only massless fermion excitations in 
the thermal bath. The same is true if we take in (|5.8| ) the limit /x ^ T as commented 
above. In our picture of chiral condensation, the Coulomb term in ( |5.8| ) represents the 
density of "residual" charges that are not yet paired to form chiral "molecules" . 

For arbitrary values of z, the solution to ( |5.6| ) has to be found numerically. First, one 
truncates the continued fraction ([5.5|) to a given order, making sure that the difference 



between successive truncations becomes smaller than a fixed numerical precision (see |22 
for details). This procedure is meaningful only if the continued fraction is convergent 
and indeed, by computing it in this way we check numerically the arguments about 
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Figure 6: The total net fermion density p(yu) (above the /i axis) and the Coulomb gas 
correction pcilj) to the massless gas (below the axis) for g'^ /tx = 9 and different values 
of the temperature around T = m. The massless result Pm=o(/^) = /^/(^r + g^) is also 
displayed for g^/Tr = 9. 



convergence discussed above. Then one has to find, also numerically, the smallest positive 
zero of Q{zo,'j,r]) = 1. For that purpose we have used a bisection method. In fact, since 
our objective is to get the function 'j{z), it is simpler to look for the zeroes in 7 fixing 



zo 



Besides, we have made use of the asymptotic expression 



to check the 



answer in the large z limit. 

Our results for pcifJ') and p{fi) are displayed in figure ^ for temperatures around 
T ~ m at fixed g'^/Ti = 9 [X'^/Att = 0.1). When z increases, or equivalently when the 
temperature decreases, pc{p) (the curves below the axis in figure H) has a similar shape as 
that given by the second term in ([5. 8]), but its magnitude increases so that the (negative) 
correction to the massless gas becomes more and more important. In fact, we see that 
for p<m, pc{p) roughly equals the massless contribution, thus yielding p{p) ^ in that 
region, whereas for large p, pc{p) becomes negligible and only the massless contribution 
remains. 

In view of this behaviour, let us define two "critical" points: /xq — "^i such that 
p{p<Po) ~ 0, p{p > po) 7^ 0, and pi such that p{p > pi) ~ Pm=o{p), where the curve 
for p{p) catches up with the massless one, i.e. pi is the chiral restoration point. Recall 
that we have no true critical behaviour here, so that po and pi can be defined only 
approximately. 

Physically, Pq ai T = would represent the energy needed to excite a fermion in the 
thermal bath. Thus the behaviour of p{p) in figure ^ near p = pq is typical of a massive 
Fermi gas: It is approximately zero until /i ~ po- For instance, notice that for a free 
gas of fermions of mass m, po{T = 0) = m and p{p, T = 0) = {p — m)6{p — m). For a 
small but nonzero temperature, that characteristic step-function behaviour is smoothed 
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because there is also thermal energy available, and this can also be observed in figure |^ 
as T increases. Remember that T ~ m is the lower temperatures we can take in our DR 
regime — see our previous comments about the problems with the T —>■ 0^ limit within 
our approach. In addition, notice that in our case, the almost exact cancellation between 
Pc'(/i</io) and Pm=o(/^^/Uo) for fow T takes place after a nontrivial numerical calculation 
of the Coulomb part, and therefore it provides a good consistency check. 

Let us now comment about the behaviour we observe for /i > /ii in figure |^. Fol- 
lowing our previous description of chiral symmetry restoration in terms of pairing of 
condensates, our results suggest that the chemical potential also increases the effective 
"coupling charge" holding the condensate pairs together. In fact, from the viewpoint of 
the Coulomb gas, we now have another way of increasing the mean potential energy with 
respect to the mean kinetic energy; namely to switch on a strong electric field \S\ ^ T. 
It seems natural then that in this regime the systems behaves as in the low temperature 
phase of Lenard, forming "molecules" . This corresponds to the /i > /^i regime in figure 
^, i.e. moderate temperatures T c::^ m and large chemical potentials fi ^ T. 

As T increases further, the effective mass of the fundamental fermion decreases with 
T, until, for very large T, the chiral symmetry is effectively restored and then p{fi) grows 
linearly with yU V/i > as given by (|5^ ). 



This type of behaviour also takes place in QCD at finite baryon density ps [2^|. There, 
at T = one has /iq — rriB — 930 MeV, the vacuum to nuclear matter transition point, 
with psifJ^o) — 0.16 fm~^ the density of nuclear matter, and at pi ^ 1300 MeV (for 
T = 0) the chiral symmetry is restored (vanishing quark condensate). At both points it 
is expected that the phase transition is of first order at T = 0. 

Another system where one observes a similar behaviour oi p{p) is the massive Schwinger 
model P^, where for large p the system is in the deconfined phase, unlike the massless 



case, where there is confinement for any p Q. 

6 Conclusions 

In this paper we have considered the massive Thirring model in 1+1 space-time dimen- 
sions at finite temperature T and chemical potential p. We have shown that in a certain 
regime, which we have denoted dimensional reduction, the statistical mechanics of this 
system is the same as that of a classical Coulomb gas in one spatial dimension, where the 
unit charge grows linearly with the temperature. 

The range of validity of the DR regime depends on the observable under consideration: 
for the pressure, the fermion density and the fermion condensate, we have seen that it 
works both for T ^ m with 5^^ > and for T>'m with g'^/ir ^ 1 (strong coupling regime). 
However, for the two-point correlators of fermion chiral operators (J±{x), it only works for 
large g, requiring in addition that the spatial distance q between the two chiral operators 
is such that q ^ (3/211. 

Thanks to the analogy with the Coulomb gas, we have been able to calculate exactly 
the pressure, fermion density, fermion condensate and two-point correlators of the MT 
model. Our results show that the chiral symmetry is restored both for high T and high 
/i in a continuous way (no phase transition). Chiral symmetry restoration takes place 
only in the strong coupling (large g^) regime. The symmetry restored phase corresponds 
to the Coulomb "condensed" phase, in which H — pairs of charges tend to pair, form- 
ing "molecules" , whereas our low T phase is the Coulomb "plasma" phase, in which the 
charges are free. The low T and high T behaviours are reversed in our case with respect to 
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the Coulomb gas because the couphng between charges grows with T. We have also seen 
that the role of the Coulomb charges is played here by the a± {x) operators. These oper- 
ators tend to pair in the high T or /i phase forming chiral invariant cr+cr_ combinations. 
These play the role of the "molecules" , whose "size" is determined by the condensation 
or screening length, which we have estimated for different values of the coupling constant 
and temperature. We find that the screening length decreases with T. In the restored 
phase, this screening length is of the same order of the inter-particle distance, which is 
proportional to the inverse of the fermion condensate. 

The case of a nonzero fermion chemical potential /i can also be related to the ID 
Coulomb gas, by noting that /i plays the role of a purely imaginary external electric 
field. For high T, the system behaves with /i as the massless case, which is consistent 
again with chiral restoration since the effective mass falls off with T. For T ~ -m, we 
see that the fermion density is very small until /io — ^ti and then starts growing (typical 
massive Fermi gas behaviour) until eventually, for large /i (/x > Hi), it reaches again the 
massless linear behaviour (chiral restoration). This is again consistent with the idea that 
the fermion excitation of the system becomes massless both at high T and at high /i. In 
the case of moderate temperatures, chiral condensates are held together at large /i by 
means of the strong external electric field. 

Among the analogies with QCD we have found, it is worth mentioning that the chiral 
condensate vanishes for both large T and /z, although in our case there is no phase 
transition. In addition, the type of behaviour of the fermion density p{n) we have obtained 
at low temperatures is also quite similar to that of the baryon number density pB in QCD, 
where po and pi are the vacuum to nuclear matter and chiral restoration critical points, 
respectively. 
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A Dimensional reduction for the two-point correla- 
tors 



In order to compare the aa correlators ([4.8|) - (^4.9|) with the corresponding // ones ( ([4.6| )- 
( [4.7|) ) in the Coulomb gas, the former must be dimensionally reduced through ( p.4|) just 
as in the case of the partition function. All calculations are performed with respect to 

Let us start with the H — correlator. The relevant integral to analyse now is: 



F2+-(A,T,L;X) = n / dr, dq,X{[Q\x, - x^] 

. _ -I J -/ u^A 



Q\xj-X) 



ejX'^/iTT 



where X = (g, r), < g < L and we have set Y = 0. The above integral is the 
counterpart of ( p.2|) . Let us consider first the limit X^/Att <^ 1. We readily realise 
that the X-dependent factor in the integrand does not contribute to leading order, since 
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^Cj = 0. Thus, following similar steps as in section p?2| , we get 
< Tca+{X)a.{0) > - < Tccr+{X)a^{0) >'^^s"" A^ m^ 



<Tca+(X)a_(0)> 



^^ lim ?7J(r?) + O — 



where /(r^) is given in ( |2.9|) and the superscript asym means now replacing the Q's inside 
the integrals as (|2.4|), but not for the function Q'^{X) outside, which would be allowed 
only for q 3> (3/2tt (see our comments in the main text). Notice once more, that one 
should be careful when taking T to arbitrary small values at the same time as expanding 
in A. Here, we even find a term m? /T"^ to leading order. Again, we will consider our 
approach valid only for T>m. 

Our numerical analysis of I{r]) shows that it vanishes much faster than r/^^ as r^ ^ oo, 
so that we see that for small A, DR (for the integrals) is justified. Notice that the above 
relative error is scale independent, so we have fixed p = m. That is not true for the 
correlator ( [4. 8] ) itself, where there is an explicit p^'^ 1'^'^ overall dependence. 

Next, we will discuss the T '^ m case. Following the same steps as for the partition 
function, we get now 



<Tca+(X)a_(0)>= 




2\ I-A^'/Stt 



ZAiriT) \p 

2 



I \ f p^ ^ 



Q\x) 



-A2/47r 



' 7T) \ 8^ A /TV 



with 



ft{\X) 



dTidT2 / dqidq2 



Q\qi 



q,n 



r)Q\q: 



'2, 1-2 



Q^(<?i - 92, n - T2)Q^{q2 -q,T2- T)Q^{qi, Ti) 



(A.l) 



X'^/iTT 



(A.2) 

where X = (g, r), g has been rescaled q — > q/P and all the Q^ functions are evaluated at 
(3=1. Notice that we have replaced p = rra in ( |A.1| ) only for the scale-independent part, 
keeping the scale dependence outside the curly brackets. 

Thus, for T 3> m, we have to analyse the behaviour of the integral ( |A.2| ) with A and X 
or, rather, the difference between /2^~(A, X) and f2~'"''^^"^{^, X) obtained by replacing the 
Q's by their asjTiiptotic values. Notice that 72*^ (A, 0) gives exactly the same contribution 
as the partition function evaluated to the same order, that is, /2^~(A, 0) = /2(A) in 
( p.8|) . In fact, in the language of Feynman diagrams, this is a disconnected contribution, 
proportional to 7], so that the factor Zq/Zmt in front ensures that the answer for this 
correlator is finite as 77 ^ cxd and we take it into account by subtracting the X = 
contribution to the correlators. The same is true for /2 ~'"*^™(A,X). Bearing this in 
mind, let us concentrate on the values of A^ close to and below Att. The reason is twofold: 
On the one hand, we have seen already that for small A the DR works. On the other 
hand, those are the values that will give us the biggest contributions to the error, since 
the integrand of the asymptotic form is not singular at any point, whereas in ( [A.2| ), there 
are singular regions where the denominator vanish, and those give larger (but finite) 
contributions to the integral the closer we approach to A^<47r. The direct numerical 
evaluation of ( |A.2| ) is a hard task. However, for T ^ m, all we need to show is that the 
difference with its asymptotic value remains bounded close to A^<47r. Therefore, let us 
estimate it by looking only at the biggest contributions, i.e, those of the regions close 
to the singularities of the integrand. Firstly, in the region xi — X2, where Xi = (gi,ri) 
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and X2 = {q2, T2), we see that the integrand simply tends to its X = value, and hence 
that contribution cancels with the partition function. On the other hand, near X2 — X 
or xi ~ 0, the integrand goes like [Q^(X)/Q^(x2 — X)Q'^{x)] ''. Then, 



ft{\X)-ft{\0)^\Q\X) 



A2/47r 



[/2(A)]^ 



(dominant contribution) (A. 3) 



We have seen in section ^ that |/2(A) - 2^'/"^^ /2\'^\<1 for A^ < in. 

As for the asymptotic contribution to the integral ( |A.2|) , we can evaluate it explicitly 
by integrating over the four separate regions on the (gi, ^2) plane according to the different 
relative signs oi qi — q2, qi — q and g2 — q- We obtain in the 77 — >• cxo limit 



f: 



H — ,asym 



(A,X)-/2+ 



,asym 



(A,0) 



2A2/27r+l c o 



A2 



-3^+A^ 



X^q/2 , i -3A2g/2 

9 



10 



(A.4) 
Therefore, from ( |A.4| ), ( |A.3D and (|A.1|) , we see that the relative error for the two-point H — 
correlator — which is a scale independent quantity — does not remain bounded at large T 
for A^ close to 4n, unlike the pressure, but it grows arbitrarily large for g ^ 1. In 
particular, taking A^<47r and g ^ 1 in the above expressions, we find 



< Tca+(X)a_(0) > - < Tca+{X)a_{0) > 
<Tc.a+(X)a_(0)> 



'1 

y2 



-~c-^exp[27rg] + -- 



2^2 



withc= (/2(A)V4-l/7r2)/4. 

Now we turn to the ++ correlator. As before, we shall start with the limit A^/47r ^ 1. 
From (14. 91) we have now 



«Tca+(X)a+(r)» 



Z!{T) 



• T-i \ -^ /27'" / „\2 CO 



^MT 



iT)\p^ 



E 



1 



^on!(n + 2)! 



m2 /T\^'/^- 



Q\q.r) 



A2/47r 



F++(A,T,X) 



2n+2 



(A.5) 



with 

F++(A,T,X) = /5^ 

where now 



2n+2 „i 
2{2n+2) T[ I ^^, 



n 







dgj Q'^{xj — X)Q'^{xj 



ejX'^/in 



Yl [Q'^{xj -Xk) 
k<j 



ejek\^/i-!T 



^3 



+ 3 
- 3 



l,...,n 
n + 1,.. 



2n + 2 



and, as before, we have set q — >• g//3, p = m, in the scale- independent part, X = (g, r) 
with < g < ?7 and the Q's in the integrand are evaluated at /5 = 1. Notice that now 
Y,j ^j 7^ 0, so that we find an additional contribution when expanding in A^/47r. We find 
for the relative error 



< Tca+(X)(T+(0) > - < Tca+(X)a+(0) > 
<Tca+(X)a+(0)> 



asym 



)2 

— lim 

At: v-'oo 



-I{v)-2I{v,X) 



Hv) 



-(£)■ 



with z = ■m'^r]/T'^, 
T] Uo 



dqi / dn log Q^(gi - g, Ti - r) + 2r] log 2 - vr?]^ + 2nq{r] - q) 
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and 






z d 
2dz 

z 

2 



.z 



/oW 



3 d 

dz 



2\2 



In the thermodynamic hmit t] —>■ oo, it is easy to check that I{ri,X) —>■ I{rj) because of 
translation invariance, and that g{z)/ f{z) — ^ z/2, so that the relative error for small A 
tends to zero (for T>m). Again, asym means the same as for the H — correlator, i.e, 
replacing the Q's by their asymptotic values only inside the integrals, but not for the 
[Q2(g,r)]^'/'" factor in (Q. 

Finally, we will consider T ^ m for this correlator. We have now 



< Tca+(X)a+(0) > 



X 




S Q'ix) 



X^/4-K 



/2++(A,X)0 - 



m 



S-AVtt 



with 



f^^{\,X) = / dTidT2 / dqidq2 
Jo Jo 



X 



Q'^iQi -q2,ri -T2) 



-1 AV47r 



Q'^iqi -q,Ti- r)Q2(g2 - g, Tg - T)Q'^{qi, Ti)Q'^{q2, T2) 



Notice that there is an important difference between this correlator and the previous 
cases we have analysed, namely, that the difference between the actual value and the 
asymptotic limit shows up already to leading order in the m/T expansion, which is now 
0{m/TY~^ /'^'^ rather than 0{1). That means that the relative error for the correlator 



/2 (A, X) starts now at C(l), and therefore is not bounded. To check that this is indeed 
the case, it is enough to take the dominant contribution to the integral by picking up the 
different poles, as we did before with the H — correlator. One easily finds now 



ft*(\X)~\Q\X) 



-"'"[m? 



(dominant contribution) 



Notice that in this correlator there is no disconnected contribution, which is consistent 
with the fact that its large T expansion begins at NLO whereas that of the partition 
function starts at 0{1). In other words, /^^(A,X) is finite as r^ ^ 00. On the other 
hand, the asymptotic contribution yields 



/2++'"'*^'"(A,X) 



2 

'A2 



2ge" 



-^^^ + 3- 



-X^q 



A2 



4- 



-WI2 

"a^ 



Therefore, one readily checks that the relative error for /2"*'(A,X) is indeed 0{1) for 
arbitrary g, so that we cannot bound it in the T ^ m limit. 

Summarising, the structure of the correlators ( [4.8|) -( ^l9D makes it possible to reduce 
them dimensionally in the small A limit but not in the large T one. As commented above, 
DR will also work for the factors Q'^{X) outside the integral in the limit q ^ /3/27r. 
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